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1. Introduction 

This article gives a brief introduction to the On-Line Encyclopedia of Integer Sequences (or 
OEIS). The OEIS is a database of nearly 90,000 sequences of integers, arranged lexicograph- 
ically. The entry for a sequence lists the initial terms (50 to 100, if available), a description, 
formulae, programs to generate the sequence, references, links to relevant web pages, and other 
information. 

2. To consult the database 

Since 1996, an electronic version [20] has been accessible via the Internet, at the URL 

ht tp : / / www . research . att . com / ~nj as / sequences / . 

If a list of numbers is entered there, the reply will display the entries for all matching sequences. 

For example, suppose you were trying to count the ways to insert parentheses into a string 
of n letters so that the parentheses are balanced and there are at least two letters inside each 
pair of parentheses. The outer pair of parentheses is to be ignored. For n = 3 and 4 there are 
respectively 3 and 11 possibilities: 

n = 3 : abc, (ab)c, a(bc); 

n = 4 : abed, (ab)cd, a(bc)d, ab(cd), (ab)(cd), (abc)d, 
a(bcd), ((ab)c)d, (a(bc))d, a((bc)d),a(b(cd)) . 

Further work shows that for n = 1, . . . , 5 the numbers are 1, 1, 3, 11, 45. Entering these into 
the web page produces nine matching sequences, but they are sorted, with the most probable 
match appearing first. Indeed, this entry tells you that these are the numbers (sequence A1003) 
arising from "Schroder's second problem", and are also known as "super-Catalan numbers". 
The reply (an abridged version is shown in Figure 1) gives 21 references, ranging from 

Schroder (1870) [17] to articles published electronically in the last few years. There is an 
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a recurrence: 



a(l) = a(2) = 1 



(n + l)a(n + 1) = 3(2n - l)a(n) - (n - 2)a(n - 1) , n > 1 



programs to produce the sequence in Maple and Mathematica, and much more. 

There's no other reference work that will carry out this kind of search. 

The Encyclopedia can also be consulted via email. There are two addresses. Sending email 
to sequences@research.att.com with a line in the body of the message saying 



will trigger the same search that the web page performs, only now the results are sent, almost 
immediately, via email. Superseeker (superseeker@research.att.com) carries out a more sophis- 
ticated analysis and tries hard to find an explanation for the sequence, even if it is not in the 
database. If the simple lookup fails, Superseeker carries out many other tests, including: 

• applying over 130 transformations to the sequence, including the binomial, Euler, Mobius, 
etc. transforms [1], and checking to see if the result is in the database 

• applying Pade approximation methods to try, for example, to express the n th term as 
a rational function of n (using the "gfun" package of Salvy and Zimmermann [16], the 
"guesss" program of Derksen [5] and the "RATE" program of Kratthentaler [8] ) 

• checking to see if changing one or two characters produces a sequence in the database. 

Since Superseeker carries out a nontrivial amount of calculation, users are asked to submit 
only one sequence per hour. 

The electronic version of the Encyclopedia had its origins in the books [18] (1973) and 
[21] (1995). Disk space is cheap, and the present incarnation (excluding illustrations) contains 
about 72 times as much data as the 1995 book. The history of the Encyclopedia is described 
in more detail in [19]. 
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3. Applications 

Most people use the Encyclopedia to identify a sequence, as illustrated above. It has been 
around long enough so that there is a good chance that your sequence will be there. If not, 
you will see a message encouraging you to submit it. 

Most of these applications are unspectacular, akin to looking up a word in a dictionary (cf. 
[2]). One encounters a sequence in the middle of a calculation, perhaps 

1 2 4 6 10 12 16 18 22 28 30 ... , 

and one wants to know quickly what it is — preferably a formula for the n-th term (in this 
case it is probably prime(n) — 1, A6093) or a recurrence. Successful applications of this type 
usually go unremarked. Some are more dramatic: there is a web page^ that lists several 
hundred articles that acknowledge help from the OEIS. One quotation will serve to illustrate 
this. Emeric Deutsch of Poytechnic University, Brooklyn, said in a recent email message: "... 
your database is invaluable. For example, for a certain sequence a n , using Maple I found the 
first 100 or so indices i for which is odd. Only the OEIS could tell me that the sequence of 
these z's is a known sequence related to the Thue-Morse sequence. Of course, this had to be 
followed by further reading and proof." 

The other main application is to find out the current status of work on a problem — for 
example, the search for Mersenne primes (see A43), the enumeration of Hadamard matrices 
(A7299), Latin squares (A315) or meanders (A5316) of order n, the latest information about 
the decimal expansion of tt (A796) or, better, its continued fraction expansion (A1203). 

Of course people trying to solve puzzles or IQ tests find the database useful. A5228 is a 
classic: 

1 3 7 12 18 26 35 45 56 69 83 ... 

There are also some less obvious applications. One is in simplifying complicated integer- 
valued expressions. You might, for example, have encountered the sum 

A / An + 1 \ / n + k \ 
2n-2k J \ k ' 

There are powerful methods for evaluating such sums [12], [13], but it doesn't take long to 
work out the first few terms: 1, 12, 240, 5376, and to look them up in the database. In this 
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case you would have been lucky. The reply suggests that this is sequence A6588, 4' 
and supplies, with references, the binomial coefficient identity you were hoping for. 

Another application is in proving inequalities. You might suspect that o~(n) < n^/n for 
n > 2, where a is the sum-of-divisors function (A203). If the initial terms of [n-y/ra] — cr(n) 
(where [] denotes the "floor" function) are submitted to the database, the reply suggests that 
this is A55682, and points you to a reference that gives a proof of your inequality 

I cannot resist mentioning sequence A57641, which gives the values of 

[H n + exp(F n ) log(H n )} - a{n) 

for n > 1, where H n is the harmonic number ^27=1 V** Lagarias [9], extending earlier work of 
Robin [15], has shown that this sequence is nonnegative if and only if the Riemann hypothesis 
holds! 

Although the database contains a number of sequences of both of the above types, I have 
not made a systematic search through reference works such as [7], [11] and it would be nice to 
get many more examples. 

The database can also be used to save space when referring to particular sequences. When 
introducing the Motzkin numbers, for example, instead of giving the definition, references and 
the first few terms it is simpler just to say ". . . the Motzkin numbers M n (sequence A1006 of 
[20])". 

One can also search the database for sequences that mention a particular name (Riemann, 
say), and there is a separate alphabetical index, useful for keeping track of all sequences on a 
certain topic — e.g. the entry for groups lists abelian (A688), primitive permutation (A19), 
transitive permutation (A2106), simple (A5180), total number (Al), . . .. 

In the past year the main look-up page has been translated into 28 languages, with the goal 
of making it easier to use throughout the world. The entries from the database still appear in 
English, but the headings in the replies and the error messages have also been translated. 

4. The database 

To be included in the database a sequence should be integer- valued, well-defined and interest- 
ing. 

The main sources are combinatorics, number theory and recreational mathematics, but 
most branches of mathematics are represented (e.g. A27623, the number of rings with n 
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elements), and there are hundreds of entries from chemistry and physics (e.g. A8253, the 
coordination sequence for diamond: the number of carbon atoms that are n bonds away from 
a particular carbon atom). 

Sequences of rational numbers are entered as a linked pair giving numerators and denomi- 
nators separately. The Bernoulli numbers B n are A27641/A27642. 

Triangular arrays of numbers are read row-by-row, so that Pascal's triangle gives A7318: 

1, 1,1, 1,2,1, 1,3,3,1, 1,4,6,4,1, ... . 

Square arrays are read by anti-diagonals, so that the Nim-addition table 

1 2 3 4 5 ... 

1 3 2 5 4 ... 

2 3 1 6 7 ... 

3 2 1 7 6 ... 

produces A3987: 

0, 1,1, 2,0,2, 3,3,3,3, 4,2,0,2,4, ... . 

Most well-defined submissions get accepted, since an open-door policy seems the best. The 
amazing coincidences of the Monstrous Moonshine investigations [4], for example, make it 
difficult to say that a particular sequence, no matter how obscure, will never be of interest. 

Sequences that are discouraged are those that depend on an arbitrary and large parameter: 
primes of the form x 2 + y 2 + 2003, say, whereas primes of the form x 2 + y 2 + 1 form a perfectly 
acceptable sequence (A79545). 

The Encyclopedia currently receives between 10, 000 and 12, 000 downloads per day. The 
rate of arrival of new sequences has remained constant at about 10, 000 per year for the past 
seven years, with roughly the same number of comments and updates. To keep this flood of 
information from getting out of control, people are asked to use a web form^ when submitting 
new sequences or comments. 

For most of its life the Encyclopedia has been maintained by the author, but in the past 
year a board of associate editors has been formed to help with the work. There is also a 
group of regular users who constantly send corrections and extensions, and help maintain the 
accuracy of the entries. 

'•http://www.research.att.com/~njas/sequences/submit.html 
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Even so, much remains to be done. There are more journals and e-print servers now than 
ever, and the trained eye sees integer sequences everywhere. I still discover articles in the 
library or on the web where authors have published sequences without sending them to the 
Encyclopedia. If you come across an integer sequence in your own work or elsewhere, please 
submit it to the Encyclopedia! 

Of course, accuracy is a major concern in maintaining the database. The entries in [18] 
and [21] were checked very thoroughly, and almost all the errors that have been discovered in 
those books were already present in the sources from which the sequences were taken. As the 
number of sequences has increased in recent years, it has become more difficult to check them 
all. However, the number of users has also increased, and a large number of the entries carry 
a comment that the sequence has been extended (or sometimes, corrected and extended) by 
someone. Contributors see a reminder that the standards are those of a mathematics reference 
work, and all submissions should be carefully checked. So, on the whole, users can be confident 
that the sequences are correct. The keywords "uned" and "obsc" indicate sequences that have 
not yet been edited or for which the definition is unclear. These serve both as warnings to 
users and to indicate places where volunteers can help. 

One of the pleasures of maintaining the database is seeing the endless flow of new sequences. 

1 will end by mentioning a few recent examples: 

Home primes (A37274), [6]: a(n) is the prime reached when you start with n, concatenate 
its prime factors, and repeat until a prime is reached (a(n) is defined to be —1 if no prime is 
ever reached, although it is conjectured that this never happens). E.g. 8 = 2x2x2^ 222 = 

2 x 3 x 37 -» 2337 = 3 x 19 x 41 -» 31941 ->...-► (after 13 steps) 3331113965338635107, 
a prime; so o(8) = 3331113965338635107: 

1 2 3 211 5 23 7 3331113965338635107 311 773 11 223 13 .... 

The EKG sequence (A64413), [10]: a(l) = 1, a(2) = 2 and, for n > 3, a(n) is the smallest 
natural number not in {a{k) : 1 < k < n — 1} with the property that gcd{a(n — 1), a(n)} > 2: 

1 2 4 6 3 9 12 8 10 5 15 18 14 7 21 24 16 20 22 11 .... 

Lacing a shoe (A78601), [14]: Number of ways to lace a shoe that has n pairs of eyelets. 
The lace must follow a Hamiltonian path through the 2n eyelets, and at least one of the 
neighbors of every eyelet must be on the other side of the shoe. 

1 3 42 1080 51840 3758400 382838400 52733721600 .... 
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A "bootstrap" sequence (A79000), [3]: a(n) is taken to be the smallest positive integer 
greater than a(n — 1) which is consistent with the assertion "n is a member of the sequence if 
and only if a{n) is odd." 

1 4 6 7 8 9 11 13 15 16 17 18 19 20 21 23 25 27 29 31 .... 
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Figure 1: Part of the reply when the sequence 1, 1, 3, 11, 45 is submitted to the On-Line 
Encyclopedia. Many references, links and comments have been omitted to save space. 



Reply from On-Line Encyclopedia 
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Greetings from the On-Line Encyclopedia of Integer Sequences ! 



Matches (up to a limit of 30) found for 113 11 45 : 



ID Number : 
Sequence : 



Comments : 



References 



Formula : 



Example : 

Maple : 
Math ' ca : 

See also: 

Keywords : 



A001003 

1, 1, 3, 11, 45,197, 903, 4279, 20793, 103049, 518859, 2646723, 
13648869,71039373,372693519,1968801519,10463578353, 
55909013009, 300159426963, 16183 62158587, 8759309660445, 
47574827600981,259215937709463,1416461675464871 

Schroeder's second problem (generalized parentheses); also called 
super-Catalan numbers or little Schroeder numbers. 

a (n) - number of ways to insert parentheses in a string of n symbols. The 
parentheses must be balanced but there is no restriction on the number 
of pairs of parentheses . The number of letters inside a pair of 
parentheses must be at least 2. Parentheses enclosing the whole string 
are ignored. 

Stanley gives several other interpretations for these numbers. 
P. J. Cameron, Some sequences of integers, Discrete Math., 75 (1989), 
89-102 . 

I. M. H. Etherington, Some problems of non-associative combinations, 
Edinburgh Math. Notes, 32 (1940), 1-6. 

S. Motzkin, Relations between hypersurface cross ratios, and a 
combinatorial formula for par titions of a polygon, for permanent 
preponderance, and for non-associative products, Bull. Amer. Math. 
Soc, 54 (1948) , 352-360. 

Schroeder, Vier combinatorische Probleme, Zeit . f. Math. Phys . , 15 
(1870) , 361-376 . 

P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; 
see page 178 . 
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Bottomley, Illustration of initial terms 



H . 

E. Pergola and R. A. Sulanke, 

Schroeder Triangles, Paths, and Parallelogram Polyominoes, J. Integer 

Sequences, 1 (1998) , #98 .1.7. 
Recurrence: a(l) = a(2) = 1; for n > 1, (n+l)a(n+l) - 

3 (2n-l) a (n) - (n-2) a (n-1) . G.f.: (1/4) * (1+x-sqrt (1-6*X+X A 2) ) . 
For n > 1, a(n) = (1/n) * sum_{k = .. n-2} 

binomial (2*n-k-2,n-l) *binomial (n-2 , k ) . 
a(3) = 3: abc, a(bc), (ab)c; a(4) = 11: abed, (ab) cd, a(bc)d, ab(cd), 

(ab) (cd) , a(bcd), a (b (cd) ) , a ( (be) d) , (abc)d, (a(bc))d, ((ab)c)d. 
tl : = (1/4) * (1+x-sqrt (l-6*x+x"2) ); series (tl, x, 40 ) ; 

Sch[ 1 ]=Sch[ 2 ]=l;Sch[ n_Integer ]:=Sch[ n ] = ( 3 ( 2n-3 ) Sch [ n-1 ] 

-(n-3)Sch[ n-2 ])/n; Array [ Sch [ # ]&, 20 ] 
See A000108 , A001190 , A001699 , A000081 for other ways to count 

parentheses. Cf. A000311 , A010683 , A065096 ■ 
nonn, easy, nice 
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